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We propose a consistent setup for the holographic dual of the strongly coupled large-iV c M = 4 super Yang- 
Mills theory plasma which undergoes the Bjorken flow relevant to the quark-gluon plasma at BNL Relativistic 
Heavy Ion Collider and CERN LHC. The dual geometry is constructed order by order in a well-defined late- 
time expansion. The transport coefficients are determined by the regularity of the geometry. We prove, for the 
first time, that the dual geometry has an apparent horizon hence an event horizon, which covers a singularity at 
the origin. Further we prove that the dual geometry is regular to all orders in the late-time expansion under an 
appropriate choice of the transport coefficients. This choice is also shown to be unique. Our model serves as a 
concrete well-defined example of a time-dependent anti-de Sitter-space/conformal-field-theory dual. 



The quark-gluon plasma produced at the Relativistic Heavy 
Ion Collider (RHIC) is an ultimate state of matter whose vis- 
cosity takes the lowest value that has ever been observed. The 
plasma will also be produced with higher temperature at the 
Large Hadron Collider (LHC). An extremely small viscosity 
indicates that the interactions among the constituent particles 
are very strong. The plasma produced by a central collision 
exhibits approximately a boost-invariant one-dimensional ex- 
pansion, namely the Bjorken flow (HI]. It is a challenging prob- 
lem to construct a theory that describes such a strongly inter- 
acting time-dependent system from the first principle. 

The Anti-de Sitter space (AdS)/Conformal Field Theory 
(CFT) correspondence provides a framework that accounts 
for a strongly coupled gauge theory in terms of a dual grav- 
itational picture in higher dimensions Finite temper- 
ature has been introduced into the framework of AdS/CFT 
in Ref. J3l and the deconfinement phase of the YM-theory 
plasma has been found to correspond to the AdS black hole. 
However, the AdS/CFT correspondence for a time-dependent 
system has not been fully established yet. 

An attempt to construct a holographic dual of the Bjorken 
flow of the AT = 4 Super Yang-Mills (SYM) plasma has been 
initiated by Janik and Peschanski J5 land has been further 
pushed forward by several groups Ja 3. H S (see also ifiolo . 
Especially, the transport coefficients including the shear vis- 
cosity have been obtained from the regularity of the dual ge- 
ometry Jil |9t] . The dual geometry has been interpreted as a 
time-dependent black hole. 



However, it is quite non-trivial even to show the very exis- 
tence of event horizon in a time-dependent geometry. As far as 
the authors know, any proof of the presence of event horizon 
on the dual time-dependent geometry has not been reported. 
Furthermore, it has been claimed in Ref. Il 111 that there is a 
logarithmic singularity at the third order of the late-time ex- 
pansion in the dual geometry which cannot be removed within 
the framework of the ten-dimensional type IIB supergravity. 

In this Letter we propose a new dual geometry which is free 
from the above mentioned problem. We prove the existence 
of event horizon, as well as the absence of naked singularity 
in the bulk to all orders of the late-time expansion. 

Bjorken flow and its holographic dual. — Hydrodynamics 
tells how the stress tensor (IV, ) of the fluid evolves in time 
provided that the equation of state (EOS) and the transport co- 
efficients are given. The EOS for the present case is given as 
the traceless condition of the stress tensor because of the con- 
formality of the gauge theory. The Bjorken flow of conformal 
fluid is simple enough to solve the hydrodynamic equation. 
The solution is given as 



T TT (r) = e (r- 4 / 3 - 2 Vq t^ + e^r^ 3 + • • • ) 



(1) 



where r is the proper time of the fluid, eo is a free parameter 
which determines the initial energy density, and ijq is propor- 
tional to the shear viscosity rj: t/(t) = eo?/o (T TT (T)/eo) 3 ^ 4 . 
We follow the second-order relativistic hydrodynamics of 



(2) 



^Present address: Asia Pacific Center for Theoretical Physics (APCTP), 
POSTECH, Pohang 790-784, Korea; nakamura@apctp.org 



conformal fluid proposed by Refs. 1121 11311 in which e ( 
given by 9l fa+ 4 ( A j.- 1 ?oT n ) w jj ere T c| anc j ^ proportional 
to the relaxation time tjj and another second-order transport 
coefficient Ai, respectively: rn(r) = (T TT (r)/eo) 1 ^ 4 , 
Ai(t) = eoA? (T T r(T)/e ) 1/2 03]. Here, the stress tensor is 
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that on the local rest frame (LRF) of the Bjorken flow and it 
is diagonal. Other components of the stress tensor, T yy and 
T x±x± where y is the rapidity and x± denote the perpendicu- 
lar directions to the collisional axis, are expressed in terms of 
T TT by virtue of the EOS and the hydrodynamic equation. 

We propose the following parametrization of the five- 
dimensional (5d) dual metric on the Eddington-Finkelstein co- 
ordinates: 



ds 2 



r 2 a dr 2 + 2 drdr + e 2b " 2c (1 + rr f dy 2 
+ r 2 e c dx 2 , . 



(2) 



Here, r is the fifth dimension with r — ► oo corresponding to 
the spatial boundary. The boundary conditions at r — > oo are 
set so that the 4d part of the metric becomes r 2 times the 4d 
Minkowski metric on the LRF of the Bjorken flow: 



ds 2 



= r 2 {-dr 2 + T 2 dy 2 + dx] 



2 drdr, (3) 



namely a — ► 1 and 6, c — ► 0. We assume that the bulk met- 
ric depends only on r and r to ensure the boost invariance 
(^-independence) fl]]. We have also assumed that the 4d part 
of the bulk metric is diagonal so that the diagonal 4d stress 
tensor is reproduced. Physically, these conditions provide in- 
formation to the gravity theory that the plasma undergoes the 
Bjorken flow. We note that there remains a gauge degree of 
freedom r — > r + /(r) in the metric (|2j. 

a, b, c (hence the 4d stress tensor) are determined from 
the vacuum Einstein equation with the negative cosmologi- 
cal constant (A = —6 in our unit), which is equivalent to (the 
bosonic part of) the equations of motion of lOd type IIB su- 
pergravity under the assumption of constant dilaton and con- 
stant Ramond-Ramond flux 112311 . Eq. (fl} suggests that the 
Einstein equation may also be solved order by order in the 
t -2 / 3 expansion. A natural trial suggested by Ref. JHt] is 
that we introduce a new spatial variable u = rr 1 ^ 3 and per- 
form a r~ 2 / 3 expansion regarding that u is independent of r: 

a(u, t) = a (u) + ai(u)T~ 2 / 3 + a 2 (ii)r~ 4 / 3 H . (This is 

similar for b and c.) We shall find that this expansion (which 
we call late-time expansion in this Letter) works consistently. 

With the boundary conditions (O, we find the following so- 
lutions to the Einstein equation: 
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where ^i, w, Ci are the integration constants which are not 
determined by the boundary conditions. We have also ob- 
tained the second-order solution, which is too lengthy to be 
shown here and will be presented in Ref. 1 1-411 - The second- 
order solution contains two more integration constants £2 and 



C2- Actually, we can show that £1 and £2 are gauge degrees of 
freedom which can be absorbed by the following coordinate 
transformation: 

u -> u - Z- 1 ^- 2 ^ - 3- 1 6t" 4/3 + 0{t- 2 ). (5) 

There is also a gauge degree of freedom of constant shift u — > 
u + £q, but we have already chosen a gauge £0 = in ©. 

Through the standard AdS/CFT dictionary lfl5tl . we find 
that the integration constants w, £1 and (2 are related to the 
energy density of the 4d fluid T TT by 



T — 

± rr 



3w 4 
16ttG 5 



T~ 4 / 3 - 2ClT~ 2 + C2T" 8/3 



(6) 



Other non-zero components T yy and T x±x± are expressed in 
terms of T TT by virtue of the Einstein equation at the vicinity 
of the boundary [141] in a consistent way with hydrodynamics. 
Namely, the hydrodynamic equation and EOS are encoded in 
the Einstein equation and the dual geometry I24I1 . Comparing 

Eqs. (O and ([T), we can read off the physical meaning of the 

(2) 

integration constants as (1 = r/Q and £2 = e o • 

Regularity of dual geometry to all orders. — The 
Kretschmann scalar computed from our metric to the 
first order is (R^ pa ) 2 = 40 + 72w 8 /u 8 + 0(t- 2 / 3 ). 
Obviously there is a singularity at the origin u = 0. We can 
also compute a component of the Riemann tensor projected 
onto an orthonormal basis 



R y 



w 



+ 0(r- 2 / 3 ), 



Ci 



4u 3 



3(3w 4 



(7) 



where = -2- x / 2 (l, 0, 0, 0, 1 + 2~ 1 r 2 a) is a fifth- 
dimensional component of the fiinfbein. The expression (0, 
which is a coordinate scalar (and a local Lorentz tensor), is 
singular at u — w except when 

Ci = l/3itf. (8) 

Hence, the spacetime regularity at u = w requires ((8). 

Similarly, we can show that the regularity of Rt^yi at u = 
w is achieved at the second order if IU4I1 



C 2 = (l + 21og2)/18u; 2 



(9) 



.-(2) 



Recalling the correspondences £i = r)o and (2 = eg ', one 
can see that the regularity of the spacetime has completely 
fixed the shear viscosity and the combination of the second- 
order transport coefficients i25ll . Let us assume that the static 
relation between the energy density T TT , the entropy density 
s and the temperature T 11711 is still valid due to the local 
thermal equilibrium which is realized at a sufficiently large 



proper time: s(r) 



-T 3 (r). Then we obtain, by us- 



ing N 2 /4n 2 = 1/8ttG 5 , that w(t)/s(t) = Si^w/Ati which 
results in the celebrated formula III 811 ??(r)/s(r) = 1/4-7T, as 
a consequence of the regularity condition (H). Similarly we 



3 



get from the second-order results that Ai(r) — t](t)tjj{t) = 
(Tt t(t) I s) , which is consistent with the results 

So far, the conditions ^ and (O are the necessary ones 
for the spacetime regularity at u = w. Actually, we find that 
these are the sufficient conditions, too. This can be shown 
by confirming the regularity of the metric and its inverse as 
well as their arbitrary-order derivatives with respect to r and r . 
(The apparent divergence of the metric at r — > oo is harmless 
as it is the case for the pure AdS.) Once these are shown to be 
regular, any coordinate invariants made of Riemann tensors 
and their covariant derivatives are also regular. The foregoing 
regularity conditions are equivalent to the regularity of a, b, c 
and their arbitrary-order derivatives in our parametrization (O 
except at the origin. 

We can generalize the above observation to all orders by 
means of induction [fl4n . The outline is the following. Let 
us use the terminology "regular" for the regularity except at 
the origin. We begin with the assumption that ak, bk, Ck and 
their arbitrary-order derivatives are regular for k < n. We 
also assume that the 1/u expansions of ak, bk, Ck around the 
boundary start at the order of 1/u or less singular order. Then, 
a component of the nth order (oc r~ 2 "/ 3 ) Einstein equation, 

(it b' n )' = terms including ak, bk, Ck (k < n) and 

their derivatives with regular coefficients, (10) 

tells us that the left-hand side is regular, where ' denotes d/du. 
We can generalize this statement to the regularity of b n and 
its arbitrary-order derivatives for < u < oo by integrating 
and/or differentiating the equation ( [Tol l. We can also prove the 
regularity of a„ and its arbitrary-order derivatives in a similar 
way, by using another component of the Einstein equation: 

(it a n )' = terms including ak, bk, Ck (k < n), b' n , and 

their derivatives with regular coefficients. (11) 

The proof for the regularity of c„ is more complicated since 
we encounter a potential singularity: 



1 



— (regular terms) . (12) 



It is necessary to prove that one can take the integration con- 
stants so that the terms denoted as "regular terms" in Eq. (fT2l 
vanishes at u = w. We find that the integration constant f„ 
in a n appears linearly in the regular terms of Eq. (fT2l hence 
it is always possible to remove the singularity by an appro- 
priate choice of £„. We also find that ( n corresponds to the 



?ith order contribution to the stress tensor e 



(») 
o 



We have ex- 



plicitly seen this to the second order: the transport coefficients 

(2) 

?7o and e are determined by requesting the regularity of c\ 
and c 2 , respectively. Since we have already shown that our 
starting assumption is valid to the second order, the regularity 
(under the appropriate choice of the transport coefficients) for 
all orders is proven by induction. 



So far we have proven that the spacetime can be made 
regular by an appropriate choice of the integration con- 
stants that corresponds to the transport coefficient . 



Now we shall prove that if we do not choose the con- 
stants £„ as above, we necessarily encounter a singular- 
ity at u = w. This is shown by computing R^-4 V a = 

u){~u if re 8 u l ar terms ^ . Any choice 
( fTZt singular at u = w necessarily renders 



7 2 T - 2/3„4( T -2/3 +l 

which makes Eq. 
R-±y± singular. 

We have proven that a unique appropriate choice of a set 
of transport coefficients makes the spacetime regular (except 
at the origin). Therefore, the unremovable singularity pointed 
out in Ref. Il 111 is absent in our setup. 

Apparent horizons. — Since we still have a physical singu- 
larity at the origin (u = 0), we shall show the presence 
of an event horizon which covers it. We show the pres- 
ence of an apparent horizon whose existence necessarily leads 
to the existence of an event horizon outside l20ll . Let us 
expand the location of the apparent horizon as uh(t) 



U + U\T 



-2/3 



U 2 T' 



-4/3 



0(t~ 2 ) and determine the coef- 
ficients order by order by solving the equation = Q(uh) = 
6 + 6it~ 2 / 3 + e 2 T~ 4 / 3 + 0(t~ 2 ), where 9 is the normal- 
ized product of expansions in the double null formalism l2lll . 
We obtain 



e = -(9/2)(i- Uo - 4 w 4 ), 

9i = -(3/t«)(l-2^i+6u 1 ) ) 

37r + 41og2-8-24w (£ 2 



e 2 = - 



3m 2 ) 



Aw 1 



(13) 
(14) 

(15) 



where the conditions for the regularity (H)) and (O have been 
imposed. The apparent horizon is located at uq = w, u\ = 
-± + ^ and u 2 = 8 ~ 3 T,~ 4 u los2 + We see the presence 
of the trapped region > inside the apparent horizon. The 
contribution of £1 and £ 2 can be absorbed by the coordinate 
transformation (O. 

The area element of the apparent horizon is given 
by Andydx^ 

m(2+7r+61og2 



where A 1 



w 



.-2/3 



24 



H - 2 

T ~ 4 /3 _|_ 0(t~ 2 ). The area element of the event 



horizon should be the entropy density per unit rapidity rs, 
that means, ts > Ah /4G '5. (At the late-time limit, the ap- 
parent horizon must coincide with the event horizon.) The 
entropy density per unit rapidity is given by ts oc 
1 — — 7p-T~ 2 / 3 + • • ■ from Eq. ([T]i and the dimensional analy- 
sis. The entropy production rate coincides with the expansion 
rate of the area of the apparent horizon at the leading order: 
Aj^drAn = (ts) _1 <9 t (ts) = r] a T~ 5 / 3 . In other words, the 
ratio of the first and the second terms in the t~ 2 / 3 expansion 
agrees between Ah and ts. 

The area element of the apparent horizon has been found to 
increase with r. Further, we can show that the difference be- 
tween the Weyl tensors exhibits a time-dependent anisotropy: 



C x y x i y = 3r] w*u- 4 T- 2 / 3 + <3(t- 4 / 3 ). This 



indicates that the geometry is not locally static under the pres- 
ence of the dissipation due to the viscosity. 
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Summary and discussion. — We have proposed a gravity 
dual of the Bjorken flow of the TV = 4 SYM fluid. The exis- 
tence of event horizon has been shown for the first time in the 
geometry dual to the Bjorken flow. We have also proven that 
the spacetime can be made regular to all orders of the late-time 
expansion by an appropriate choice of the integration con- 
stants corresponding to the viscosity, the relaxation time, and 
other transport coefficients. Furthermore, this choice is not 
only sufficient for the regularity but also necessary, namely, 
the requirement of the regularity completely fixes these coeffi- 
cients. The area element of the apparent horizon has also been 
computed. The increasing rate of the area element coincides 
with that of the entropy density per unit rapidity at the leading 
order. Our spacetime exhibits a time-dependent anisotropy in 
the Weyl tensors, indicating that the geometry is not locally 
static. 

It is interesting to consider what classifies the proper- 
ties determined around the boundary and those determined 
around the horizon. The EOS and the hydrodynamic equa- 
tion have been obtained at the vicinity of the boundary while 
the transport coefficients have been obtained from the regu- 
larity around the horizon. The traceless property (EOS) and 
the conservation of the stress tensor (hydrodynamic equation) 
hold whether or not the (local) thermal equilibrium is achieved 
in the YM-theory side. However, the concept of the transport 
coefficients makes sense only when the notion of fluid is valid. 
This tempts us to relate the notion of (local) thermal equilib- 
rium with the regularity (or the presence) of the horizon. It is 
also interesting to see how the method to determine the trans- 
port coefficients from the regularity is related to the Kubo's 
linear response theory (Kubo formula) and other holographic 
computations. We hope that these points will be clarified in 
the future. 

We have discussed the time evolution of the entropy density 
from the viewpoint of the dual geometry. It would also be in- 
teresting to study a thermodynamic formulation of dynamical 
black holes in asymptotically AdS^ geometries. 

The present paper provides a consistent setup for the holo- 
graphic dual of the Bjorken flow of the Af = 4 SYM 
plasma. Our model also provides a concrete well-defined ex- 
ample of time-dependent AdS/CFT. We hope that the present 
work sheds some light on the dynamical nature of the time- 
dependent plasma. 

Note added: When the present work was at the final stage, 
we received a paper lEill which partly covers a related sub- 
ject. The first-order solution presented in Ref. 12211 corre- 
sponds to the gauge choice of £1 =0 in our first-order so- 
lution. Our proposal in the present paper has been invented 
independently. However, we were motivated by Ref. 1E2I1 to 
examine the gauge degree of freedom and the regularity of the 
higher-order geometry. 
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